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Abstract 

An accurate and comprehensive room temperature rotation-vibration transi- 
tion line list for phosphine ( '^PHa) is computed using a newly refined potential 
energy surface and a previously constructed ah initio electric dipole moment 
surface. Energy levels, Einstein A coefficients and transition intensities are 
computed using these surfaces and a variational approach to the nuclear motion 
problem as implemented in the program TROVE. A full ro- vibrational spectriim 
is computed. The resulting line list, which is appropriate for temperatures up 
to 300 K, consists of a total of 13.7 million transitions between 2.56 million en- 
ergy levels covering the wavenumber range to 8000 cm"-'^. Several of the band 
centres are shifted to better match experimental transition frequencies. The 
line list is compared to the most recent HITRAN database and other labora- 
torial sources. Transition wavelengths and intensities are generally found to be 
in good agreement with the existing experimental data, with particularly close 
results for the rotational spectrum. An analysis of the comparison between the 
theoretical data created and the existing experimental data is performed, and 
suggestions for future improvements and assignments to the HITRAN database 
are made. 



Preprint submitted to Elsevier 



February 11, 2013 



1. Introduction 

Phosphine (PH3) is an extremely toxic, rigid molecule found in the lower 
troposphere of the Earth and is expected to be of great importance for the 
spectroscopic studies of giant-planets, particularly for probing the deeper layers 
of the atmospheres [1], and cool stars. Phosphine has for a long time known 
to be abundant in Jupiter [H El H], where it is partially responsible for the 
continuum opacity in the 5 /xm region of the atmosphere [5^ and where it is a 
marker for vertical convection zones. It has also been found on Saturn [HH], 
and is expected to be present in extrasolar gas giants. 

Further spectroscopic investigations of phosphine in the atmospheres of as- 
trophysical bodies requires an extrapolation of the spectra over a wide range 
of temperatures, pressures and path lengths. The ability to accurately identify 
and interpret phosphine features in astrophysical spectra is dependent on the 
existence of a complete description of the phosphine spectra, for all relevant 
temperatures. 

The current JPL [6] and HITRAN [7] databases encapsulate decades of ac- 
curate laboratorial measurements of phosphine spectra at room temperature. 
However, between the two databases, only a total of 20827 lines are described, 
all within the — 3600 cm~^ region, with J < 23. These databases are designed 
to be valid for temperatures below 300 K. Other sources (see Table [l] below) 
mostly overlap with the HITRAN and JPL databases, and only add a little 
more to these compilations. 

The quantity and variety of spectral information required to correctly char- 
acterise hot astrophysical bodies is often beyond what can be expected to be 
delivered experimentally. The ExoMol project [S] (see www.exomol.com) aims 
to provide a solution to this problem by performing accurate quantum mechan- 
ical calculations, with the aim of providing appropriate spectroscopic data on 
all relevant molecules for the analysis and characterisation of cool stellar and 
exoplanet atmospheres. This is, however, a challenging alternative to measur- 
ing spectra, requiring highly accurate potential energy surfaces (PES), dipole 
moment surfaces (DMS) and an ability to generate precise wave functions and 
eigenvalues for the nuclear Schrodinger equation and the associated transition 
intensities. In practice, at least for the PESs, it is also necessary to incorporate 
a certain amount of experimental refinement. 

Our aim is to develop complete line lists which could be used to accurately 
model the atmospheres of stars and Jupiter-like planets, for which a high tem- 
perature (> 1500 K) line list for PH3 is of special importance. In order to 
establish and test the production of such a line list, which will consist of ap- 
proximately 1 billion transitions, it is first pertinent to perform a simulation of 
the room-temperature intensities for this molecule, which can be tested against 
existing reliable experimental sources. In the present paper a cool (300 K) line 
list is presented. 

There have been several previous attempts at ab initio studies of the behavior 
and characteristics of phosphine spectra. Wang et al [3] computed a PES using 
coupled cluster CCSD(T) theory and Dunning's correlation consistent cc-pVTZ 
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basis sets, followed by a vibrational analysis based on second-order perturba- 
tion theory. They determined several spectroscopic constants and fundamentals 
largely within 4 cm~^ of the observed results. Zheng et al [10] calculated an 
ab initio three-dimensional P-H stretching DMS using density functional theory 
and found absolute band intensities agreed with observed ones within a factor 
of 2. Yurchenko et al [TT] presented a calculated, albeit empirically refined, 
PES and, using a variational method, calculated the rotational energy levels in 
the vibrational ground state of PH3 for J < 80 [l2]. Subsequently, Yurchenko 
et al |13j computed an entirely ab initio six-dimensional dipole moment sur- 
face (DMS) at the CCSD(T)/aug-cc-pVTZ level of theory for the electronic 
ground state of PH3 . This is the DMS that is used in the present work, and is 
described in detail below. It gives calculated transition moments within an av- 
erage of 22.5% of the experimental data. The same paper also presented a PES 
determined by empirically refining an existing ab initio surface. Ovsyannikov et 
al [13] complemented the variationally calculated PES and DMS of Yurchenko 
et al [13] by calculating a PES at the CCSD(T) level using aug-cc-pV(Q-hd)Z 
and aug-cc-pVQZ basis sets for P and H, respectively, and then presenting a 
list of computed vibrational transition moments for the electronic ground state 
of PH3 [15]. The PES used here is a refined version of the PES presented by 
Ovsyannikov et al [T3j, and is discussed further below. Recently, Nikitin et 
al produced a new PES [TB] and used it to calculate several vibrational en- 
ergy levels with only a marginal deviation from experiment. Nikitin et al then 
successfully modelled the lower three polyads of phosphine using a variational 
approach jl7^. The resulting line list is discussed and compared to the present 
work in Section 4.6. The phosphine line list presented here can be considered 
the analogue of the ammonia line lists computed by Yurchenko et al |18| , using 
TROVE. 

2. Background 

Phosphine is a well behaved symmetric top molecule belonging to the group- 
15 hydrides. The phosphorous atom is positioned on the axis of symmetry, 
perpendicular to the plane containing the equidistant three hydrogen atoms. 
As is common with molecules having Csv (M) symmetry [T^ , there is a splitting 
of rotational levels with K — 3n (where n > 1) in non degenerate vibrational 
states [in]. Phosphorous has only one non-synthetic, stable isotope, '^^P, so only 
^^PHs is considered in the present work. 

Phosphine has near degenerate vi symmetric and 1/3 antisymmetric stretch- 
ing modes [5] with frequencies of about twice that of the asymmetric bending 
mode, V4. It is therefore natural to consider its spectrum in terms of polyads, 
which will be described in section 3.1. There is also a strong Coriolis interaction 
between the z/2 and 1^4 bending bands which causes large distortion effects in 
observations 21] . Recorded phosphine spectra show no inversion splitting. The 
tunnelling effect found in ammonia is expected to occur in phosphine but so far 
attempts to detect it have failed due to its considerably higher potential energy 
barrier (12300 cm^^) between the two symmetry-related minima [22] [23]. 
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In the absence of the inversion sphtting, phosphine is characterized by the 
C3v{M) molecular group symmetry, i.e. each eigenfunction transforms according 
to one of the irreducible representations Ai, A2 and E, the latter of which is 
two-fold degenerate. 

To fully describe the internal rovibrational motion of phosphine one needs 
the following minimal set of quantum numbers |24| : 

JT-i , "2 , "3 , "4 , ^3 J -^4, Tvib , J, K, Trot , Ttot , (1) 

where L3 = |Z3|,L4 = 1^41, i = \^iK = \k\. Here the vibrational quantum 
numbers ni (symmetric stretch), 712 (symmetric bend), (asymmetric stretch) 
and 714 (asymmetric bend) correspond to excitations of, respectively, the i^i, 1^21 
z/3 and modes. The doubly degenerate modes and require additional 
quantum numbers L3 = 1^31 and L4 = |^4| describing the projections of the 
corresponding angular momenta (see, for example, Ref. [H]). The vibrational 
quantum number L = |/|, characterizes the coupling of Z3 and Z4. Fj-ot, Tvib, 
and Ftot are, respectively, the symmetry species of the rotational, vibrational 
and total internal wave- functions in the molecular symmetry group C3v(M), 
spanning Ai, A2 and E. J is the total angular momentum and K — \k\,k = 
—J, . . . , J is the projection of the total angular momentum on the molecule fixed 
axis z. 

Out of these twelve quantum numbers, only two are rigorously conserved 
quantum numbers (J and Ftot)- This quantum number description is similar to 
that proposed by Down et al for ammonia. The rigorous dipole selection 
rules include transitions between states with Ftot satisfying Ai ^ A2 and E ^ 
E and A J = 0, ±1 only. Strong transitions obey the dipole selection rules 
A J = ±1 and /S.K = 0. There are also weakly allowed transitions which 
occur when, due to centrifugal distortion, the C31, geometrical symmetry of the 
molecule is broken and K is no longer a good quantum number |25| . 

Only approximate selection rules can be associated with the vibrational 
quantum numbers ni, n2, 713, ri4. 

Table [T] presents an overview of the experimental spectra recorded for phos- 
phine. This spans research starting from 1951 to the most recent body of work, 
published in 2009. The naming convention for each individual source consist 
of two digits for the date and first two letters of up to four authors' surnames, 
e.g. S. Yurchenko, M. Carvajal, P. Jensen, H. Lin, J. Zheng & W. Thiel 2005 
becomes 05YuCaJeLi. 
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Tabic 1: Experimental sources of phosphine spectra. 



Source No of lines Range(cm ^) Uncertainties(cin ^) Intensities 
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3. Method 



The production of the phosphine hne list presented here relies on the com- 
putational program TROVE (Theoretical RO Vibrational Energies) [H], which 
uses a general variational approach to calculate the ro-vibrational energies for 
small semi-rigid polyatomic molecules of arbitrary structure. It employs a nu- 
merical finite basis representation. The nuclear kinetic energy operator is nu- 
merically constructed through a recursive scheme using a Taylor series expansion 
in terms of the molecule's internal coordinates, which requires no analytical pre- 
derivation, making its creation self-contained. This process is an extension of 
the XY3 program [45 . Although this procedure allows TROVE to simulate the 
nuclear motion for an arbitrary molecule, it does not calculate an exact kinetic 
energy operator. As such, the accuracy of the rotation-vibration energies de- 
pends on the expansion order of the kinetic energy operator and its associated 
level of convergence. Similarly, the potential energy operator is expressed as a 
Taylor-type expansion, and its numerical derivation can lead to an accumulation 
of round-off errors. For a more detailed description of the TROVE functionality, 
including its treatment of the kinetic and potential energy operators, the reader 
is referred elsewhere [H]. 

3.1. Convergence 

To construct the ro-vibrational basis set, the contraction scheme described 
perviously [JHl 071 HH] was followed. This method employs 1-D basis functions 
0i(?f ) (« = 1 ■ • . 6), where , and ^| are the linearized versions of the three 
stretching coordinates Ar^, i — 1,2,3 (Ar^ = rpHi — Te), and ^5, and 
are the linearized versions of the three bending coordinates Aa^, i = 1,2,3 
associated with the variation of the bond angles 023, ais, and ai2, respectively. 
Here fe and ae are, respectively, the equilibrium values for the bond length and 
the inter bond angle facing the bond length. 

The primitive basis functions for each mode Ui, i — 1 . . . 6, are generated 
by solving the corresponding Schrodinger equation with the Numerov-Cooley 
method [321 110] • In order to control the size of the basis sets at all contrac- 
tion steps the polyad-truncation scheme [33] is used, based on the PH3 polyad 
number defined by 

P = 2(si + S2 + S3) + 61 + 62 + &3 < -Pmax (2) 

where Si and bi are the quantum numbers associated with the basis functions, 
4>si and </>(,. , for the stretching modes and the bending modes, respectively. The 
basis set is formed only by those basis functions 4>s, 4>b for which P < Pmax- 

The primitive basis functions obtained at this step are then used to solve 
two 3-D Schrodinger equations for each class of modes, stretching and bending, 
separately. The corresponding basis sets are formed as direct products of the 
primitive basis functions satisfying 2(si + S2 + S3) ^ Pmax (stretching) and 
61 + 62 + &3 < Pmax (bending) in agreement with Eq. ([2]), so the maximal vi- 
brational excitation is taken to be -Pmax/2 and Pmax for each of the stretching 
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and bending modes, respectively. The resulting two sets of eigenfunctions 
and are then classified using C3^(M) symmetry by analysing their trans- 

formational properties and assigned with the normal mode quantum numbers 
rii , rig ^ (stretching modes) and n2 , n^^ (bending modes) for future ease of line 
assignment and comparison to the experimental spectroscopic data. To this 
end a correlation between the primitive and normal mode quantum numbers 
is needed to be constructed. At the initial J = stage, this is a straightfor- 
ward process due to the relatively small number of solutions and simple polyad 
structure of PH3 energies: only approximately 700 sets of vibrational quantum 
numbers have to be translated for the full 2.56 million energy levels (J < 31) to 
be represented by both sets of quantum numbers. 

The reassignment from local mode quantum numbers i^i [i = 1 ... 6) to 
normal mode quantum numbers Ui [i — 1...4), L3 — [l^l and L4 = {l^l was 
performed by applying the following rules : 

S1+S2+S3 = ni+ns, 

bi + b2 + b3 = n2+ni, 

h = -"-3,-"3 + 2,...,n3 - 2,n3, 

I4 = — ri4, — 714 + 2, 714 — 2, 714. 

It was also assumed that the symmetric modes and 1^2 have lower energies 
than their asymmetric equivalents 1/3 and 1^4, respectively, and also that the 
vibrational energy grows when I3 or I4 increase. 

At the next contraction step, a 6-D vibrational basis is formed as a direct 
product of the stretching and bending functions and <&JiJ°"^, contracted using 
the (normal mode) version of the polyad number truncation given by Eq. ([2|: 

P = 2(711 + n^) + 772 + 774 < -Pmax- (3) 

The 6-D functions (|)bend ^^^^ then symmetrized by the standard re- 

duction technique [T^. At the next stage the vibrational ( J = 0) Schrodinger 
equation is solved variationally by diagonalizing the (J — 0) Hamiltonian matrix 
on this 6-D vibrational basis. Using the Pmax polyad restriction from ([3| 1455, 
1125, and 2571 basis functions are obtained for the Ai, A2 and E symmetries, 
respectively. 

At the last contraction step, the resulting eigenfunctions of the 6-D J = 
problem are used to form our final ro-vibrational basis functions in the so- 
called J = representation |46) through a direct product with the symmetrized 
rigid rotor wave functions \J,K,t), where r = 0, 1 indicates the rotational 
parity (—1)^ as defined in Ref. [45]. In the J = representation building the 
Hamiltonian matrix is straightforward as the vibrational part of the matrix is 
diagonal. 

The ro-vibrational eigenfunctions obtained can be naturally assigned using 
a scheme based on the largest contribution of the basis set expansion. The 
chosen quantum numbers are the vibrational normal mode quantum numbers 
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ni, 712, 71,3 ^ and n^'^, the rotational quantum numbers J, K, the vibrational sym- 
metry Trot, the vibrational symmetry Fvib, and the total symmetry Ftot, where 
Trot/vib/tot — Ai,A2,E. The projection of the total vibrational angular mo- 
mentum L is also included to reduce ambiguity in the description of the energy 
levels, with the following rules: L must be a multiple of 3 if Fvib = ^1,^2 and 
cannot be a multiple of 3 if Fvib = E. By convention, the lower energy value was 
assigned to IL3 — L4I. However, at higher energies, TROVE does not necessarily 
assign unique quantum labels to every state. In this case we have chosen to 
assign all states to IL3 — L4I. 

The accuracy with which high ro-vibrational states can be computed depends 
heavily on the size of the Hamiltonian matrix. To ensure that the calculation 
of the energy levels was sufficiently complete, convergence tests were performed 
to facilitate the choice of the polyad number (Pmax), ranging from P = 4 to 
P = 18. A selection of the results from the convergence tests is given in table 

m 



Table 2: J = energy eigenvalue convergence with polyad number, Pmax, used to generate 
the basis set, where P(n) = (Pmax = 18) - (Pmax = 



Band 


■Pmax 18 


P(16) 


P(14) 


P(12) 


P(10) 


P(8) 




0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


l>2 


992.732 


-0.010 


-0.033 


-0.090 


-0.243 


-0.714 


2U2 


1967.427 


-0.133 


-0.393 


-0.939 


-2.197 


-5.321 


2ua 


2222.257 


-0.009 


-0.033 


-0.095 


-0.283 


-0.897 


V\ 


2324.077 


-0.003 


-0.011 


-0.036 


-0.128 


-0.436 


3U2 


2911.866 


-1.325 


-3.472 


-7.230 


-14.335 


-29.377 


V2 + 2V4, 


3204.920 


-0.115 


-0.361 


-0.941 


-2.467 


-7.220 


V\ + 1^2 


3322.008 


-0.038 


-0.132 


-0.384 


-1.135 


-3.820 




3341.959 


-0.028 


-0.098 


-0.295 


-0.931 


-3.434 


1^3 + 


3439.304 


-0.009 


-0.034 


-0.117 


-0.437 


-1.843 


^-V2 


3803.487 


-9.750 


-22.284 


-39.911 


-66.534 


-108.521 


2V2 + 21^4 


4159.298 


-1.025 


-2.839 


-6.309 


-13.381 


-28.564 


V\ + 2V2 


4303.413 


-0.352 


-1.074 


-2.740 


-6.742 


-16.622 


V2 + 3!^4 


4318.966 


-0.257 


-0.786 


-1.981 


-4.906 


-12.643 


+ + i'4 


4410.959 


-0.131 


-0.436 


-1.196 


-3.221 


-9.642 


4!/4 


4422.841 


-0.098 


-0.325 


-0.919 


-2.707 


-7.461 



Ideally, the polyad number Pmax would be chosen to give full convergence, 
which, as can be seen from table [2] is not achieved even for the highest polyad 
configuration, Pmax = 18. However, although the calculations should improve 
as the polyad number increases, since the 'spectroscopic' PES (see Section 3.2) 
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used was generated with Pmax = 14, the energy levels move away from the 
observed values when a higher polyad configuration is used. The PES used in 
this work relies on a self consistency of parameters, and as such it is only an 
effective potential for a specific set of parameters, so Pmax = 14 was chosen. 

Convergence tests were performed to test the expansions of the kinetic and 
potential energy operator and, although accuracy did improve considerably as 
the expansions grew, so did the computational cost involved. The convergence 
test presented in Table [2] was done using expansions of the kinetic energy op- 
erator and potential energy function both to fourth order, which accounts for 
the significant discrepancies between the computed band origins and the cor- 
responding experimental values. In the present work, the expansions of the 
kinetic energy operator and potential energy function are truncated at the G*'' 
and S"* order, respectively, which leads to much more accurate results, as can 
be seen below. Higher order expansions than these would have made the project 
computationally prohibitive. 

3.2. Potential Energy Surface 

The PES of phosphine used here is a refinement of the ah initio (CCSD(T) / aug- 
cc-pV(Q+d)Z) PES [13], done by performing a least square fit to available 
experimental ro-vibrational energies with J = 0,1,2,4 and 10. These refine- 
ments followed the fitting procedure introduced in Ref. [47]. The refinement 
of the PES is represented as a correction /S.V to the ah initio PES, Vq. The 
ro-vibrational J = 0, 1,2,4, 10 eigenfunctions of the 'unperturbed' Hamiltonian 
Hq — T -'r are used as basis functions when iteratively solving the set of 
Schodinger equations for the ro-vibrational behaviour of the molecule to min- 
imise the associated functional in the least-squares fitting. With the final refined 
PES, a total root-mean-squares (rms) error of 0.03 cm^^ was obtained for the 
fitting energy set. Table [3] shows the change in the band centre values between 
the pre- and post-refinement PES. 

As already mentioned, the kinetic energy operator is not calculated exactly 
and the basis set is artificially limited. Thus the refined PES must be considered 
an effective PES; it only gives the accurate results presented in Table [3] when 
used with TROVE and the parameters described above. The refined PES used 
here is given in the supplementary material to this paper in the form of a Fortran 
program. 

The PES is a major source of error in the line list creation process, and as 
such constant updates and refinements are necessary. For example, although 
the current PES is appropriate for creating an accurate room temperature line 
list with frequencies up to 8000 cm~^ and J < 31, it would not necessarily 
remain so for a line list with parameters exceeding this. For example, the PES 
fails to accurately describe the potential of phosphine at very high excitations, 
showing some artificial minima in the region of dissociation. Consequently, this 
PES is not suitable for use in dynamical simulations without adjustment. 
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3.3. Dipole moment surface and transition intensities 

An existing six-dimensional ab initio electric dipole moment (CCSD(T) / aug- 
cc-pVTZ) |13j is used to obtain the Einstein coefficients and transition inten- 
sities. This was calculated on a large grid of 10080 molecular geometries. The 
procedure used to compute absolute intensities mirrored that used for ammo- 
nia, and is described in detail in Yurchenko et al [18| . A number of vibrational 
transitional moments were calculated to help characterize the quality of the 
DMS. 

Table [3] compares the empirical vibrational transition moment with values 
computed by TROVE, both in the present work and in the previous studies 
[T3] . It is clear that the theoretical band intensities are in good agreement with 
experiment and, although the DMS used here is unchanged from that used in 
Ref. [T3j, the use of our more accurate PES (and hence wavefunctions) has 
meant that our new results have generally reduced the error in the calculations 
of the transition moments, from an average deviation of 22.5% to 10%. 



Table 3: Calculated band centres and their respective transition moments (deviation from ex- 
perimental values shown as a percentage). Uncertainties of the experimental (Obs) transition 
moments are given in parentheses (in units of the last digit quoted) where available. 



Band Centres (cm ^) Transitions Moments (D^) 



Band 


Obs 


Calc '06 13| 


Present 


Obs 


Calc '06 (%) [13] 


Present (%) 




0.000 


0.000 


0.000 


0.57395(30) 


0.583(1.6%) 


0.585(2.0%) 


1^2 


992.130 


992.500 


992.152 


0.08251(5) 


0.085(2.5%) 


0.084(2.2%) 


1/4 


1118.310 


1117.870 


1118.322 


0.08626(5) 


0.087(0.3%) 


0.089(3.7%) 


2U2 


1972.550 


1972.820 


1972.590 


0.00299(5) 


0.003(9.7%) 


0.004(21.4%) 


V2 + Vi 


2108.150 


2107.170 


2108.169 


0.01102(6) 


0.009(15.6%) 


0.014(25.3%) 


2u4 


2226.830 


2227.860 


2226.835 


0.0176(2) 


0.006(68.8%) 


0.018(1.5%) 


2u4 


2234.930 


2234.570 


2234.940 


0.0176(2) 


0.002(91.5%) 


0.013(23.5%) 


VI 


2321.120 


2322.040 


2321.142 


0.0690 


0.073(6.1%) 


0.072(4.3%) 


Ki 


2326.870 


2329.180 


2326.888 


0.130 


0.139(6.8%) 


0.138(6.1%) 



3.4. Empirical adjustment of the vibrational band centres 

Although using a 'spectroscopic' PES improves the values of the energy 
levels, the TROVE calculations do not completely reproduce the observed tran- 
sition frequencies. To correct this, an empirical approach was adopted, where 
an artificial frequency shift is added to calculated band origins as given by the 
J = energies; this procedure is described as an empirical basis set correction 
scheme (EBSC) [TH]- It leads to a rotational energy structure in much better 
agreement with experimental results for the remaining J values of the band. 
As a test, the subsequent J = 1,2 and 4 values were compared to those from 
experiments which, together with J = were also the values of J used to refine 
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the PES. While using this method, care must be taken to only select reliable 
experimental data, since data of limited accuracy can be a source of error which 
TROVE cannot compensate for. 

Only 11 band centres were manipulated, and several iterative shifts were 
attempted. The final improvements to the standard deviation, cr, of the whole 
band from experiment are displayed on table |4] The original average a value 
for the set of bands was 0.031 cm"-'^ which was reduced to 0.02 cm~^ after this 
final adjustment. 



Table 4: Observed band origins (Obs) and standard deviation with which TROVE reproduces 
the terms within each band, a, before and after replacement of the band origins. Observed 
data is from HITRAN2008 ^. 



Band 


Obs 


Original a 


New a 




992.135 


0.030 


0.020 


1/4 


1118.307 


0.012 


0.005 


2V2 


1972.571 


0.011 


0.007 


V2 + 


2108.152 


0.045 


0.034 




2226.835 


0.029 


0.010 




2234.920 


0.032 


0.014 


v\ 


2321.121 


0.028 


0.012 


vz 


2326.8667 


0.023 


0.013 




2940.767 


0.063 


0.046 


1^2 + 2l/'4 


3214.936 


0.035 


0.024 


V'i + 1^4 


3440.259 


0.037 


0.034 
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4. Results 



For the room temperature line list presented here the following thresholds 
were selected. The chosen ranges of energy eigenvalues are 4000 cm"-'^ for the 
highest lower energy and 12000 cm^^ for the highest upper energy. These choices 
allow for a range of transitions of to 8000 cm^^ at temperatures up to 300 K. 
The lower energy threshold of 4000 cm^^ defines the highest J that had to be 
taken into account, as there were many more energy values for each symmetry 
and each J than those within the selected range. By J = 31 only 0.0025% of all 
energy eigenvalues (a total of 8 levels, across all symmetries) were below 4000 
cm~^ and by J = 32 there were none. Figure [l] shows the variation in number 
of useful (within range) eigenvalues as J increases. 

Summing over all the energy levels calculated here gives a partition function 
equal to 3249.5 at 296 K. This is 0.028% higher than the latest published HI- 
TRAN value of 3248.6 [5T]- 



250 




Figure 1: Total number of energy eigenvalues under 4000 cm ^ for J = 0, 1, 33. The blue, 
green and rod lines represent, respectively Ai, A2 and E symmetry eigenvalues. 

For the purposes of comparison with the existing experimental data, only 
transitions above a minimum intensity (> 10~^ cm/mol, or 1.7x 10~^^ cm/molecule) 
were included. This is approximately one order of magnitude more sensitive than 
the weakest lines in the HITRAN 2008 database [7| . With this intensity cut-off, 
5 086 749 transitions were selected from the total of 13 725 5400 computed lines. 

It is possible to use this line list to generate synthetic spectra as a function 
of temperature. However, such spectra will become increasingly incomplete as 
the temperature goes above 300 K. 
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Figure |2] shows an overview of our simulation compared to both JPL [5^155] 
and HITRAN [7] databases. HITRAN lacks any pure rotational transitions and, 
although it captures most lines stronger than 5 x 10~^^ cm/molecule, it is very 
incomplete below this value. 
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1000 2000 3000 4000 5000 6000 7000 8000 

Wavenumber (cm"*) 



Figure 2: Comparison between simulated absorption spectrum created by TROVE and those 
from the JPL [6] and HITRAN fTl databases, as a function of frequency. 

All the information on the 2.56 million rotation-vibration energy levels and 
the associated list of the 13.7 million transitions can be found online at www. exomol. com, 
in the ExoMol format [54]. Additionally, the temperature-dependent molecular 
absorption cross sections [55] for phosphine is also available online. A HITRAN 
format line list is included as supplementary material to this paper. An inten- 
sity cut-off of 10~^^ cm/molecule is used for the rotational part of the spectrum 
(0 — 500 cm"-'^) and 10~^* cm/molecule for the remainder of the spectrum. Ad- 
ditionally, the refined potential energy surface used here is provided as a fortran 
program. 

Data for the rotational spectrum was obtained from the JPL database [S2], 
although about half of the lines in the JPL database can also be found in the 
CDMS data [ST. The documentation from JPL was released in 1991 and cites as 
its sources Helms and Gordy, [27^, Chu and Oka ^25j, Davies, Newman, Wofsy, 
and Klemperer [26] j Krupnov, Melnikov, and Skvortsov (31], Pickett, Poynter, 
and Cohen [29] and Belov et al. [2^. Further information about these papers 
can be found in Table [T] The JPL data set contains 728 transition lines in the 
region - 201 cm~^ with a maximum J quantum number of 27. 

An algorithm created with MATLAB matched theoretical lines produced by 
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TROVE with equivalent ones in the JPL dataset. These matches deviated from 
those in JPL with a rms deviation of 0.055 cm~^. Figure [s] shows close-ups 



h 1 1 1- 



26.66 26.67 26.68 26.69 26.7 26.71 26.72 26.73 26.74 
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113.8 114 114.2 114.4 114.6 114.8 115 115.2 115.4 115.6 115.8 164 164.5 165 165.5 166 166.5 167 167.5 168 168.5 169 

Wavenumber (cm"') 



Figure 3: JPL (lower) versus TROVE (upper). The matched TROVE transitions are high- 
hghted in red. 



of the rotational spectrum, where a hne by hne comparison can be made with 
the JPL data mirroring TROVE's theoretical predictions, and those transitions 
that could be matched highlighted in red. Some of the K=3n transition doublets 
are very close so their intensities are combined in the JPL database, but the 
equivalent TROVE lines are left resolved. This difference in data handling leads 
to an apparent factor of two disagreement in the intensities. 

The most recent HITRAN data [7] for phosphine contains 20 099 transitions 
in the region 770 - 3601 cm~^. Of these, only 16 008 are assigned and of these 
2011 have no upper vibrational quantum numbers, possibly due to vibrational 
mixing. There seems to be no consistent description of symmetry and many 
transitions appear either mislabelled or misassigned. 

The energy levels of the assigned lines were compared to those calculated 
in TROVE and this information was used to match lines between experiment 
and theory. Due to the ambiguity (between Ai and A2) in many of the energy 
levels with A symmetries, only the E symmetry transitions were matched with 
sufRcient confidence for an accuracy analysis, although an attempt to find a 
match for every labelled HITRAN line was made. Assigned E symmetry transi- 
tions correspond to 7838 of the total HITRAN lines. 13 of these transitions had 
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ambiguous or incorrect K values, so only 7825 E transitions were matched. This 
sample can be expected to give representative rms for each band. These were 
found to have a rms deviation from their experimental equivalent of 0.23 cm~^. 
The frequency wavenumbers of a few transitions involving energy levels which 
are defined by a single transition in the HITRAN database disagree significantly 
with that of those calculated with TROVE. These skew the rms deviation and 
its true value is expected to be significantly lower. In fact, when the top 1% 
worst matches are removed from the comparison, the rms deviation lowers to 
0.19 cm-^ 

The first region (Polyad number — 1) is located between 770 - 1372 cm~^ 
and is dominated by the fundamental bending bands vi and 1/4 and the "hot" 
band 2vi - v^. TROVE's rms deviation from HITRAN in this region is 0.23 
cm~^. Its rms is 0.37 cm~^ for the V2 band, 0.11 cm~^ for the 2^4 band and 
0.11 cm~^ for the 2v2 — vi band. HITRAN references Brown, Sams and Kleiner 
[55] as the source for this region. Figure |4] shows close-ups of a selection of 
representative sub-regions within the P = 1 region. 
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Figure 4: HITRAN (lower) versus TROVE (upper), with matched transitions highlighted in 
red, for the region P = 1. 



The second region is located between 1885 - 2478 cm^^ and is dominated 
by the fundamental stretching bands v\ and v^^ the 1vi and 2v/^ fundamental 
overtones and the combination band vi + v/^^. TROVE's rms deviation from 
HITRAN in this region is 0.20 cm-\ with 0.28 cm^^ for the band, 0.22 
cm^^ for z/3, 0.05 cm"-*^ for 2v2^ 0.15 cm^^ for 2vt, and 0.08 cm~^ for V2 -I- v^. 
HITRAN references Tarrago et al [T] as the source for this region. Figure [5] 
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shows close-ups of a selection of representative sub-regions within the P = 2 
region. 
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Figure 5: HITRAN (lower) versus TROVE (upper), with matched transitions highlighted in 
red, for the region P = 2. 



The third region is located between 2733 ~ 3600 cm~^ and is dominated 
by the4f2 — 1^2 hot band, the 31^2 overtone and the h'l + 1^2, 1^1 + 1^4, 1^2 + t'Sj 
2^3 + ^"4, 1^2 + 2i^4 and 21^2 + i^i combination bands. TROVE's rms deviation from 
HITRAN in this region is 0.27 cm^\ with 0.46 cm^^ for 4i/2 - 1^2, 0.37 cm^i 
for 3^'2, 0.50 cm"-'^ for i/i +1^2, 0.21 cm~^ for i^i +1^4, 0.23 cm^^ for 1/2 + 1^3, 0.19 
cm-i for j/3 -I- 0.11 cm-i for 1^2 + 2i^4 and 0.28 cm"! for 2i/2 + 1^4. HITRAN 
references Butler et al [ID] as the source for this region. Figure [6] shows close-ups 
of a selection of representative sub-regions within the P — 3 region. 

Nikitin et al [16] computed 55 223 transitions in the range 700 - 3500 cm~^ 
(the HITRAN range [7]). This work does not aim for completeness and even 
within the HITRAN regions appears to have significant omissions (e.g. maxi- 
mum J = 20). The work presented here improves on their root mean squares 
deviation from observed band centres of 1.4 cm~^. 

Experimental data (1768 lines) from Wang et al [37] in the 1950-2480 cm^^ 
and 3280-3580 cm~^ regions was compared to the theoretical lines created by 
TROVE. The results presented here deviated from those of Wang et al's with an 
overall rms deviation of 0.17 cm^^, for those lines that could be matched with 
confidence. The rms deviation for the first region was 0.11 cm"^, compared to 
0.20 cm~^ for the equivalent HITRAN region and 0.21 cni^^ for the second. 
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Figure 6: HITRAN (lower) versus TROVE (upper), with matched transitions highlighted in 
red, for the region P = 3. 



compared to 0.27 cm"! for HITRAN. 

Additionally, our intensities agree significantly better with the Wang data 
than with HITRAN's data, as can be seen by the intensity plots in Figure [7) 
suggesting that Wang et al's intensities should be used in a future release of 
HITRAN. 




-lo' ' ' ' ' ' ' ' ' ' ' 51 ' ' ' ' ' ' ' ' ' ' ' 

-8-6-4 -2 2 4 6 S 10 12 1 2 3 4 5 6 7 8 9 10 11 12 

log(Intobs) 

Figure 7: Plot of log{Idhs/ Icalc) versus log{Idhs) for the HITRAN data 7 (left) and the 
Wang data [37) . 
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5. Discussion and Conclusion 

In this work a phosphine line list is produced with 13.7 million transitions 
between 2.56 million energy levels for rovibrational states up to Jmax = 31 
and energies up to 8000 cm~^. This work replicates very well the observed 
phosphine spectra at room temperature, with a maximum rms deviation from 
JPL of 0.055 cm~^ for the rotational spectrum and of 0.23 cm~^ from HITRAN. 
This is a valid line list for any phosphine analysis below 300 K, however our 
ultimate aim is to construct a hot line list capable of simulating observed spectra 
from astronomical bodies at higher temperatures, e.g. 2000 K. This line list 
will complement the one already available for ammonia [46j . Work on this is 
currently in progress. 

The tunnelling effect present in the ammonia molecule is predicted to be 
found in phosphine 22 but, due to its much higher barrier (12 300 cm~^), yet 
to be observed. The value of splitting in various vibrational states as well as 
the intensity of the inversion-rotation, and inversion-rovibrational lines can be 
computed by adapting the procedure given here to work with D3;i(M) symmetry. 
Predictions for transitions which can be used to resolve the doublet splitting will 
be presented in future work. 
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